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Abstract 

£T) [ Scaling functions, F + (u/uij) and F_(lu/uj c ) for (p > <ftc and 4> < 4*ci 

O'S respectively, are derived from an equation for the complex conductiv- 

. ity of binary conductor-insulator composites. It is shown that the 

QVy \ real and imaginary parts of F± display most properties required for 



the percolation scaling functions. One difference is that, for uj/uj c < 1, 
5RF_(w/a; c ) has an w-dependence of (1 + t)/t and not lo 2 as previously 
predicted, but never conclusively observed. Experimental results on 
a Graphite-Boron Nitride system are given which are in reasonable 
agreement with the cj( 1+ *)/* behaviour for $tF^. Anomalies in the real 
Q . dielectric constant just above cp c are also discussed. 

> : 

^ : 1 Introduction 

.5?. 

The ac and dc conductivities of resistor and resistor-capacitor (R-C) net- 
works and continuum conductor-insulator composites have been extensively 
studied for many years. In systems where there is a very sharp change (metal- 
insulator transition or MIT) in the dc conductivity at a critical volume frac- 
tion or percolation threshold denoted by <p c , the most successful model, for 
both the dc and ac properties, has proved to be percolation theory. Early 
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work concentrated on the dc properties, but since it was realized that 
the percolation threshold is a critical point, and that the percolation equa- 
tions could be arrived at from a scaling relation, several papers, which are 
referenced and discussed in a previous paper 0, reporting experimental re- 
sults on the ac conductivity have appeared. Review articles, containing the 
theory and some experimental results, on the complex ac conductivity and 
other properties of binary metal-insulator systems include ||-[]nj. 

In another paper, extensive dc conductivity and low frequency dielectric 
constant results on systems based on Graphite (G) and hexagonal Boron Ni- 
tride (BN), which have what are probably the cleanest and sharpest dc MITs 
yet observed in a continuum system, have been reported O, O]. These G- 
BN results were found to obey the percolation equations, as a function of 
volume fraction in [|11|, |TI|. In [[7] it was shown that the experimental disper- 
sion results, for samples with various volume fractions of G, can be scaled 
onto two curves which are consistent with previously measured percolation 



parameters |TT|, |12|1 ; one curve refers to the real conductivity above and the 
other to the imaginary conductivity (real dielectric constant) below the criti- 
cal volume. Unfortunately, the parameters uf that had to be used to achieve 
this scaling were found to be different from those expected from the scaling 
models (R-C lattice || and anomalous diffusion no matter whether the 
uj+ and u~ were calculated using the accepted universal parameters or the 



already reported dc conductivity and low frequency dielectric results [0, 112 
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However, there is agreement between the critical exponents, characterising 
the frequency dependence of the real conductivity when > C (where de- 
notes the conducting volume fraction), the imaginary conductivity for < C 
and the exponents found from the previously reported dc and low frequency 



ac results QTT], [12 



One measurable scaling function that did not agree with the previous 
power law predictions ||-[[nj was that for the loss component in the insu- 
lating region. As these results were somewhat controversial they were not 
discussed in R. In the meantime further measurements have been made 



on one of the G — BN systems |Tl| , |12|j and other systems, using a recently 
acquired dielectric spectrometer system which is able to measure far smaller 
dielectric loss parameters as well as loss and phase angles. Some of these 
results are presented in this paper. 

The present paper introduces scaling functions which can depend on com- 
plex conductivities. They closely fit results of the medium conductance o~Mi 
in particular the frequency (u) dependence of the first order real part for 
> C and the first order imaginary part for < C as was shown in [^]. 
Most of the scaling power laws given in |H-||10| are obtained, the range of 
parameters over which these functions can be expected to generate accurate 
scaling functions is derived. It is also shown that, while the second order 
terms of the scaling functions have the exponents t/(s + 1) for « C and 
ijj jijj c ^> 1 as given in pJ-jjTO], the second order exponents for low frequencies 
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and > or < 1 differ from those of [§[|-||10|. New experimental results for 
ctm(u), 0) are presented which agree reasonably well with the exponents pre- 
dicted by the new scaling functions, provided that the complex conductivities 
of the dielectric components of the continuum systems are taken into account. 
The measurements for o~m (0 < 0c) in the G-BN systems are definitely not 
in accord with the u 2 prediction given in [p|-|10|. 



A feature of the new experimental results is that, where measurable, the 
real dielectric constant continues to increase with for > C and certainly 



does not decrease according to (0 — C )~ S as given in ||-[1(J. However this 
increase is qualitatively consistent with the expressions introduced in this 
paper. The effect is more clearly observed in carbon black-polyethylene 



compounds |T3[] and in 3D-systems where various conducting powders are 



distributed on the surfaces of large insulating grains [15 

2 Theory 

The equation 

(i-0)(^-4*) , _ n m 

gives a phenomenological relationship between ac, crj and o~m which are the 
conductivities of the conducting and insulating component and the mixture 



of the two components, respectively |7|, |T^]. Results obtained from an earlier 
version of Eq.(|l|) are reviewed in jTj| and references therein, where the 
dc conductivities of some two phase systems are successfully modelled for 



s = t. The conducting volume fraction ranges between and 1 with 
= characterising the pure insulator substance (<tm = 07) and 0=1 
the pure conductor substance (cr M = 077). The critical volume fraction, or 
percolation threshold, is denoted by C , where a transition from an essentially 
insulating to an essentially conducting medium takes place. We use the 
notation A = (1 — C )/0 C . For s = t = 1 the equation is equivalent to the 



Bruggeman symmetric media equation fig , \Vq . The equation yields the two 
limits 

S 

Wc\ 0° : °M = (Tjjj ^-T7 < 0c (2) 

(0c - <t>) 8 

I I n (<t>-<t>cY , ^ , 

Wl\ -> : M = CTc-yq TT7 > 0c (3) 

(1 - 0c)' 

which characterise the exponents s and t. Note that Eqs.(Q) and (§) are the 
normalised percolation equations. For ac measurements, || || have given 
equations for the case where ac is real and 07 = — icueoe/, which characterizes 
a lossless dielectric. However, we note that all three quantities 07, 077 and 
(7m can in principle be complex numbers in Eq. ([[]). A solution for om can 



be obtained after rewriting Eq.(|I]) using the variable z = a 



l/t 

M J 



VIZ. 



Az l+a _ + jyl/t _ ^ _ ^ _ ^a/t _ (^a^l/t = q (4) 

with a = t/s. We note in passing that Eq.(|J) has explicit solutions for 
a = 1, 2 and 3, while numerical solutions are easily obtained for larger integer 
values. Our interest is now focussed on the question as to what extent the 
solution for 07/ can be used to obtain valid scaling functions. 



The scaling conditions, which are based on those given in P, |9|, |I0|, read 
ou = ctc (0C ~ 0) V „(x„), 0<0 C (5) 



c 



° M = ° C (\ _ A\t F +( X+ ^ > 0c (6) 



where the scaling functions F±(x±) depend on the scaling parameters 

x - = — 72 277772 = < 0c 7 

°c C - S+ ^c 



with 



and 



with 



- °C (0c - <PY +t 



<?i (1 ~ 0c) s+t .UJ 
x + = — 72 — 272772 = ^ > 8 

0"C (0 - 0c) S+ * W+ 



+ _ CTC (0~0c) g+t 

Wc e o e / (l-0c) s+t ' 



The expressions involving cj± assume specifically a purely real ac and imag- 
inary 07. To ensure that curves drawn for F± fall on top of each other for 
different <f) c , the normalisation employed in all the equations used in this 



paper differs somewhat from the one used in || || |TTJ. Using the variable 
u = Ft! 1 an equation is found for u by the substitution z = uoq{^) c — 4>)/4> c 
in Eq.(£|). It reads for <p < <f> c 

A 1+a , a (0 C -0)(1 - 0-0 c ) 1/s l/s n (a s 

Au + u — u — xl — xl = 0. (9) 

0c 
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In a similar way, the substitution z = uaj ((ft — C )/(1 — (ftc) leads to an 
equation for F+ 1 (again denoted by u) for (ft > (ft c 

Au 1+a - Au« - u ^-^-f-^) x r/s _ j,. = Q (1Q) 

C (1 - <Pc) 

If the term linear in u of Eqs.(|9|) and fllPf ) could be neglected, the scaling 
functions F± would manifestly depend only on the respective variables x±. 
It is due to this term, that scaling is invalidated to a certain degree by the 
solution for <jm of Eq.(jl]). The range and extent to which this is the case is 
discussed below. An interesting aspect of Eqs.(|) and (pi]) is exact scaling 
at (ft = <p c and <ft = 1 — <p c . Whether or not exact scaling for <fi = 1 — <p c is 
merely a coincidence can only be revealed by appropriate experiments. We 
stress that all results obtained in this section are independent of whether the 
conductances and hence the scaling functions F± are genuinely complex or 
real. 

Exact solutions of Eqs.(|9]) and (|TUp can be read off at the limit points of 
the concentration. At <\> = 0, it is u a = x_ s , i.e. F_ = x_ and at = 1, it is 
u — 1, i.e. F + = 1. From these solutions, Eqs.(|2|) and (H) are obtained from 
Eqs.(5) and (6), respectively. In fact, the respective solutions are valid to 
high accuracy for > and (ft < 1 as long as \x±\ 1 or for (ft very close to 
either (ft c or 1 — (ft c . Correction terms are given below. 

For the opposite limit of the scaling parameters, i.e. \x±\ ^> 1, we obtain 
for the leading term at (ft m (ft c the solution Au 1+a = x 1 ^ from both equations, 
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(||) and (|K]). This translates into 

F± = (x±)^A~^t (11) 

which gives 

vai = ^t(— at = C . (12) 

Note that Eq.flTjp conveniently lends itself for complex values of 07 and uq', 
in particular, if 07 is purely imaginary and oc real one obtains 

= ~— n: h— l 3+t sin ( 9 / a t rr )» ( 13 ) 

J^m = +4- — ' +| cos L ■ a 14 

A~ °c 2(s + t) 

From Eq.(|ll|) it follows that the slope of the real and imaginary part of 
l°g(-F±) is t/(t+s) when plotted against log(x±) for \x±\ ^> 1. Equations (|T3|) 
and (|TJ) show that, with 07 = —iue ei, the frequency dependence of both 
real and imaginary <jm is . This dispersion law is given in f?||-||10||. Ex- 
perimental results validating this power law are found in J7L I2|,|EQ]-[23]. The 



loss angle 5 = arctan[7rt/(2(s + t))] implied by Eqs. ( JIBLpp is also given in 

It is of physical interest to determine the correction terms of next order 
in Eqs.(||) and Note that Eq.(§) yields a purely imaginary result for om 
if a i is imaginary. However, a loss term should emerge for uo > when > 0. 
This is obtained by expanding the solution of Eq.flj^) for <fi > 0. One finds for 
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1 , 1 

t / J. / Gl 

XL <Ji 



the scaling function 

F - = x -- s i x -^— = -( 1 + ( s + t ^ ~ S T2—^ — ) ( 15 ) 

which can be used for real or complex 07 or ac- Combining Eq.(jI5[) with 
Eq.(|5|) one obtains for small 0/0 c , as expected, both an enhanced dielectric 
loss term + s<p/(f) c ) and a composite loss term. Taking specifically 07 

purely imaginary and ac real, the composite loss term reads explicitly up to 
terms linear in 



07 



, r 1 sm^J 

3?ctm = sac-^ — — 1 • (16) 

Al^-I* + 2A\^\~ cos(^) + 1 

An important consequence of Eq.(^) is the small frequency behaviour of the 
loss term (recall 07 = —iueoei), which implies 

9^0" m ~ w ^ , (17) 
which differs from the a; 2 -behaviour predicted by the expansions used for F_ 
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1C]. We note that these expansions assume analytic behaviour for om 
around to — which is in contrast to our findings; also we obtain a loss term 
that vanishes for <fi — > 0, which is not the case for the expressions in [j^j-QTO] . 
The following section presents experimental results which appear to confirm 
the power law expressed by Eq.fll7D. 

By similar means we obtain the first order correction term in the vicinity 
of < 1 which reads 

f+=1+f =1+f -1=4- (18) 

C (1 - <Pc) a + xl C (1 -(pc) A + 



Note that this term implies not only a correction to the real part of <tm in 
Eq.(|3|) but also a switching on of an imaginary part for complex 07. For 
purely imaginary 07 this is 

1-6 I^H'sraff) 
StoM = -t " *. <ro ac \ 23 r (19) 

6 c (l-6 c ) A 2 + 2A|^|"cos(f ) + 1^1" 

1 ac 1 v 2s ' 1 ere 1 

which implies in this limit, for 07 = — iueoei, that 9<tm ~ cj 1 ^. 

So far, we have concentrated on regions where scaling is obeyed by the 
solution of Eq.(|l|) either exactly or to high accuracy, that is for < 6 < 6 C 
and 6 C < 6 < 1, if |x±| 1, and for 6 m 6 C , if \x±\ ^> 1. There is an 
intermediate region \x±\ ~ 1, where the linear term in u of Eqs.(^|) and 
([10]) does invalidate the sole dependence of F± on x± except for 6 = 6 C or 
6 = 1 — 6 C . In fact, it can be shown that, as long as the inequality 

e e 7 

is obeyed, the linear term of Eqs.(|9]) and ( |TUD is immaterial and scaling 



^ <Pc . , 

— < with ujq 

UJq 1 - 20 c 



(20) 



prevails. As a consequence, the leading behaviour of F±, for x± 3> 1, is 
governed by the power law x l ± t+sS} only up to the frequency which obeys the 



inequality fl20D , for larger frequencies F± becomes a linear function of x±. 
Note, however, that the right hand side of Eq.(|20D depends on 6 C in such a 
way that scaling is expected to be invalidated only for small values of 6 C and 
sufficiently large values of uj. In turn, for 6 C > 1/3 the right hand side of ( p0|) 
is larger than unity, and for 6 C —>■ 1/2 no bound on uj prevails. (Note that 
6 C = 1/3 is the Bruggeman value for spheres in three dimensions and 6 C = 1/2 
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for discs in two dimensions ||T6|, |I7|.) As a consequence, there should be no 
discernible deviations from scaling for C > 1/3. To what extent these results 
are physically valid can only be assessed by experiment. No experiments in 
this region appear to exist and the situation is complicated by the fact that 
oc and oj depend on uj when u becomes sufficiently large. 

In Fig.l we illustrate the behaviour of the real and imaginary parts of 
F±(x±) for s = 1, t = 2 and = <p c = 0.16, that is for the situation where 
scaling holds exactly. As discussed above, deviations are marginal when is 
near to C and become noticeable only when the inequality (p0|) is appreciably 
invalidated. Note the equal slopes for large x± of all four curves in accordance 
with Eqs.(|l^,|14]). Also, since t > s, the imaginary parts are larger than the 
real parts; for s = t all four curves would coincide asymptotically. For s > t 
the real parts would be larger than the imaginary parts but no such system 
has been observed or predicted. 

We interpret the dependence on aj/ac of the percolation loss term as 
predicted by Eqs.( |l5| - |T7| ) as follows: Consider a three dimensional lattice, 
with ideal capacitors on nearly all the bonds but with a small number of 
randomly distributed resistors embedded in the capacitor matrix, either as 
isolated components or as small clusters, and choose Cuo < 1/R. If a voltage 
is applied to two opposite faces of the lattice, the displacement current in 
the lattice is determined almost entirely by the capacitors. Therefore, as 
the current in the resistors is evoked by a "fixed current source", which is 
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determined by the value of l/(Cu), the dissipation is proportional to I 2 R, 
i.e. the larger R the more power is dissipated. This argument can be ex- 
tended to continuum systems, which qualitatively explains the dependence 
on oijoc in Eq . (|i~5"f ) . The presence of the exponent t, which determines 
the rate of increase of the conductivity beyond (f> c , would therefore appear 
to have a role in the formation of the conducting clusters, which determine 
the dissipation below C . However, below C , the behaviour of the complex 
effective conductivity is dominated by the imaginary component, which is 
primarily determined by the interconnectivity of the insulating medium, or 
more specifically its ability of keeping the clusters of the conducting medium 
disconnected below C , which in turn is characterised by the exponent s. 

We note that the treatment given in ||S||-||10|| has no t— or s— dependence 
for x_ < 1 or x+ < 1, respectively. However, our Eqs. (|i5|) and fll8|) which are 
based upon Eqs.(^) and (]I0D indicate that the exponent t, which characterizes 
the formation of the conducting backbone, continues to play a role for <fi < <ft c ; 
similarly, the real dielectric constant is not independent of s for <p > <p c . 

3 Experimental Method 

The percolation system |7], |ll|, |12[ which best exhibits the (l + t)/t behaviour 
is a lightly poured powder of 55%g Graphite 45% Boron Nitride which is com- 
pressed, expelling air, in a capacitive cell through the percolation threshold. 
As the percolation threshold is at 0.124 (volume fraction of G), the insulator 
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at and around this point consists of 11.4%BN and 88.6% dry air, which obvi- 
ously has a 3?e close to one and a very low term. The dissipation in more 
compacted systems, such as compressed pellets of G-BN |7j or polyethylene- 
carbon black |]T8[| , would appear to be dominated by the dissipation in the 
dielectric component. 

As the constructions of the cells and experimental procedures for the G- 
BN powder have been adequately described for the dc and low frequency ac 
measurements in [JTlJ, 0] as well as the dispersion measurements (<ta/(u;) for 
fixed (j>) in JF| they will not be repeated here. Some experimental results 
are presented, which are obtained by tumbling various conducting powders 
with larger wax coated insulating grains. The conductor coated grains are 
then compressed into discs ]19| . 



The new measurements presented here are results of the real and imagi- 
nary parts of the conductivity between 10 _1 and 3.10 6 Hz, obtained using a 
newly acquired Novocontrol Dielectric Spectrometer. This instrument is able 
to measure far smaller loss components of the dielectric or insulating phase 
(equivalent to a resistor of 10 18 f2 at 10 _1 Hz and 10 8 Q at 10 5 Hz in parallel 
with a perfect capacitor) and had a better resolution of loss or phase angle (a 
maximum of tan 5 of > 10 3 and a minimum of < 10~ 3 can be measured) than 
the instruments used in R or measurements of a similar nature pOfl-pT 
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4 Results and Discussion 



The experimental results for Jie and JJctm for the 55%G-45%BN powder as 
function of frequency between 1CT 1 and 3.10 6 Hz are given in Figs. 2 and 3. 
The dispersion results for the three highest lying curves in Fig. 3 are conduct- 
ing samples (0 > <f> c ) with dispersion free conductivities at low frequencies 
[j§|l-||10|. The dielectric constant of these three samples (upper curves in Fig. 2) 
show a strong dispersion, which should go from s/(s+t) = 0.47/ (0.47+4.8) ~ 
0.09 at cf) c to 1/s ~ 2 near to <p = 1. The observed values range from 0.1 
to 0.33. The dielectric constant for the conducting samples (<f) > <f) c ) will be 
further discussed below. 

The sample marked by plus signs in Fig. 3, (0 — C = —0.0004 ± 0.001) is 
actually metallic, as its conductivity breaks away from a constant slope at 
sufficiently low frequencies, as predicted for metallic samples in this paper 
and f|-[[nj. Note too that the larger exponent 0.15 for the higher frequencies 
shown by this sample in Fig. 2 also indicates that it is not an insulator (x_ < 
1). When compared to the results in Fig.l in J7|, this observation illustrates 
the necessity of making measurements at lower frequencies than previously 
done in any experiments of this nature. 

In Fig. 2, the exponent of 9fte drops form -0.10 to -0.02 for the next four 
samples. While, from the results given in J7j, a dispersion exponent of -0.10 
could indicate a sample still in the x_ > 1 region, the fact that the sum of 
the absolute values of the exponents for the dispersion in 3?e and ^Rxjm are 
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greater than one, precludes this possibility f?||-fl0|| . This, combined with the 
ever decreasing magnitude of the 3?e exponents, going down to 0.02 from the 
lowest lying curve, allows us to conclude that these four samples lie in the 
dielectric or x_ < 1 region. Note that the lowest measured values for s/(s + t) 
recorded in the literature f/jpU-Pfl are 0.07 for 50%G-50%BN and 0.10 for 
55%G-45%BN [|. 

Based on the arguments given above, the four lowest lying plots of the 
conductivity against frequency in Fig. 3 are the conductivities for samples in 
the insulating or dielectric state (x_ < 1). There is a very slight upward 
curvature of the results between 10 and 3.10 6 Hz, but the mean slopes in 
this region are 1.06 and 1.09. The slope decreases below 3Hz. According to 
Eq. flT5| ) the conductivity is made up of both a dielectric and percolation loss 
term, the relative contributions of which must now be examined. 

In order to evaluate Eqs. QTB|,[TB|) , (JTE|) and ([0]), an expression for 07 is 
required. Unfortunately the loss term in the Boron Nitride- Air system is too 
small to measure directly on the dielectric spectrometer, which also meant 
that the percolation parameters for this system could not be measured di- 
rectly. Therefore, the complex conductivity of Boron Nitride-Air mixture had 
to be calculated, using effective media theories, from measurements made on 
a compressed disc with a porosity of 0.19. When plotted against the fre- 
quency, in the range 10 _1 to 3.10 6 Hz, 9fte for this disc was found to be virtu- 
ally constant, and that the ?R,it m term could be fitted to 3?o"m = tT dc + Du> 9 
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[ f2~5|j with o"dc being sensitive to how dry the BN in the disc was (and therefore 
also how dry the BN in the G-BN powder was). By pumping on the discs 
long enough the dc component became unmeasurable (i.e. < 10~ 18 (f2m) -1 
at 10 _1 Hz). D was independent of a^ c and had a value of 9.0 x 10~ 16 . Un- 
fortunately, the powders could not be pumped to lower values of <Td c as this 
caused them to collapse. 

Although the system is anisotropic 0, if measurements are made in the 
axial (compression) and radial directions, it is probably still valid to use 
the Hashim-Strickman (H-S) upper and lower bounds [f2f| for measurements 
made in the radial direction only. Therefore, dtaM and 3<tm for "bulk" BN 
where determined using the upper bound, as the BN grains are obviously 
in contact at low volume fractions, which makes the system closer to one 
where the BN surrounds the air [27]. The bulk parameters are 9fte = 4.1 and 
D = 1.2 x 10 -15 . As the critical volume fraction for the 55%G-45%BN system 
is 0.124, the upper limit for dilute systems was taken to be = 0.02. At this 
volume fraction of G the volume fraction of BN was 0.0164 and of air 0.9636. 
Therefore the parameters determining <jm for a system with 0.017 volume 
fraction of BN and 0.983 for air were calculated, again using the formula for 
the H-S upper bound, which gave a 3?e = 1.03 and a D = 1.41.10 with 
Ode selected to fit the experimental results. The calculated Sicr^ at (f> = 0.02 
is then used to evaluate the enhanced dielectric loss term, the percolation 
contribution and the combination of these, using Eqs.(|5|) and (|T5D. These 
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are all shown in Fig.3. The computed value of 9?<tm for the Air-BN insulator 
is not shown but lies below the enhanced curve by a factor of s0/0 c - The 
other values s = 0.47, t = 4.8 and ac = 3126(ftan) -1 , used to calculate these 
curves, have all been determined by dc percolation measurements JTl], [12| and 
Ode was chosen to match the curvature at low frequencies. The mean slopes 
of 1.05 for the squares and 1.09 for the crosses are lower than the value of 
( 1+4.8) /4. 8 = 1.21, as is predicted by Eq . (|17D . The discrepancy is explained 
by the comparatively large values of (f) in the experimental results; in fact 
numerical solutions of Eq.(f|) show that the cjVCW) behaviour is flattening 
out for increasing (ft as it has to attain the cu*/( s+i ) behaviour for <p — > <p c . 

The only data that claims an uo 2 dependence for ctm, 4> < 4> c , is that of 
Benguigui However there are problems with his data in a mixture of 

iron balls and glass beads. His sample contains only 10 5 particles (about 
10 8 in the present experiments not even counting the air volume) and only 
30 particles between the capacitor plates (nearly 300 plus air in the present 
experiments). Although Benguigui mentions the dc conductivity of the glass, 
it would appear that neither this nor the dielectric loss term in the glass is 
taken into account. The authors are also at a loss to explain how the very 
low frequency dispersion observed for dte(u) in an insulating 20.0% iron balls 
sample (Fig.2 in |ZT]) can give a ^te{u)uj ~ Scr^u;) which varies as uj 2 when 
plotted against the frequency (Fig. 8 in Therefore we do not regard 

these experiments as definitive. 
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In Fig. 2 the dielectric results for the conducting samples are terminated 
at 1 Khz as for lower frequencies tan S exceeds 10 3 and the dielectric spec- 
trometer gives spurious results. However, the results at 1 Khz (a frequency 
commonly used in low frequency experiments to measure die) clearly show 
that 3fte continues to increase with <p below and above <fi c . This is in sharp 



disagreement with the predictions given in [j8[1 — 1|10|| , where it is claimed that 
!fte should decrease as (0 — 4> c )~ s for <fi > <p c . 

The smooth behaviour of 3?e as a function of <fi passing and extending 
beyond </> c has also recently been observed in carbon black-polyethylene com- 
posites |0| and a number of systems where various fine conducting powders 
are impregnated onto the surface of almost spherical insulator grains before 



the coated grains are compressed into a three dimensional continuum [|1S 
Therefore there is now strong experimental evidence that the second order 
term for <fi > (j) c given in P|-|10| is in disagreement with experimental evi- 
dence. Their second order percolation term also fails to vanish for cf) — > 1 in 
contrast to our result in Eq.([T8|). 

However, Eqs.(|l]) and show a 3?e that continues to increase above C . 
Unfortunately, the agreement with the experimental results is qualitative, at 
best, if the parameters obtained from dc experiments are substituted into 
Eq.(|]), as is shown in Figs. 4 and 5. 

Figure 4 shows the experimental 9fte results at 1 Khz and 1 Mhz, plotted 
against 0, for the G-BN powder. All theoretical curves, calculated from Eq. 
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(H) are for t = 4.8 |lTfl . The solid curves are for oc = 3.1 • 10 3 (fim) _1 
[0, pi]] , 3?e = 1.18 (calculated C value), C = 0.124 and s = 0.96. The upper 
solid curve is for 1 Khz and the lower one for 1 Mhz. The dotted curves 
use the same parameters except that s = 0.6. As the experimental results 
all lay above the theoretical curves, we display the dashed curve where oc 
has been changed to 3.1 • 10 5 (fim) _1 . From this figure and the behaviour 
of the theoretical plots it is apparent that the experimental results can be 
better fitted, if some or all of the above parameters are varied to get the best 
fit. In this case it would not be necessary to change <p c from its dc value of 
0.124 ± 0.001. As the dc conductivity changes by nine orders of magnitude 
between 0.120 and 0.130, one cannot argue that the incorrect C has been 
identified. 

Of the six powders that were coated onto the wax coated insulating grains, 
as previously described, only the Nickel powder did not show an increasing 9fte 
above the dc value for C , this was because the conductivity at and above <p c 
was too high to measure 9fte [19|. This combined with the results for carbon 



black-polyethylene [18], both function of 0, and C ( cLS cL function of 

temperature) for fixed <p, show that this could well be the usual behaviour for 
continuum systems. The reason that this has not been previously reported 
is probably due to the limitations of tan 5 in the instruments previously 
used, and the fact that for t w 2 and large oc /(^eo e ) ratios the solution of 
Eq.(|4]) is sharply peaked just above C . To observe this difference accurate 
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measurements of both the dc and the apparent ac value of C (delayed peak 
in die) would have had to be made. This has never been done except for the 
G-BN system, the coated grain system [l9j and to some extent the carbon 



black - polyethylene system Hl8| , Pqj . In all the coated grain systems sharp 
changes in the dc curves of &m give unambiguous C values . With a poor 
conductor (magnetite) 3?e increases sharply near C and then smoothly up to 
30 c El. 

Figure 5 shows the experimental results for 3?e for a system of insulating 
grains coated with NbC ||I9|| . All theoretical curves, calculated from Eq.(|4]) 
are for t = 4.8. The solid curves are for ac = 7000(f2m) -1 , 9fte = 7.45, C = 
0.065 and s = 0.8. The upper solid curve is for 1 Khz and the lower one 
1 Mhz. The lower dotted curves are for the same parameters but s = 0.4, 
which is the measured dc value. The low value for ac could be due to the 
resistance of the system being largely determined by contacts between the 
extremely hard and angular NbC grains. All the above parameters, except 
for s = 0.8, are close to those obtained from the dc conductivity fit using 
Eq. ([!]). Except for the fact that the drop in e for larger values of is not 
observed the results could be called qualitatively correct. 

The curves selected do not show that the theoretical curves widen con- 
siderably as t is increased, i.e. the conductivity exponent plays a large role 
in determining 9fte. An explanation may be that as the conductivity of a 
percolating system above C increases more slowly with for high t values, 
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this means that a smaller fraction of the conducting component is on the 
backbone of a system with a higher t value than a system with a lower t. 
This "off the backbone" conducting material then creates nearly conducting 
links between different sections of the backbone, which are broken by the 
insulating component. It would then be the effect of these inter dead end 
capacitances that continues to increase 3?e above C , until these capacitances 
are shorted out by an ever more conductive backbone. 

5 Summary and Conclusions 

In this paper we have shown analytically that Eq.(|T|) in conjunction with 
Eqs.([5]) and (0) gives complex scaling functions for continuum percolation 
type systems with the following results: 

i) 3lF + (x + ) has zero slope for x + < 1 (first order term), and a slope of 
t/(s + t) for x + > 1, 

ii) ^F^(x_) has slope unity for x_ < 1 (first order term), and a slope of 
t/(s+t) for x_ > I. 

iii) A slope of t/(s + 1) for $sF± and ?R.F± for x + and x_ > 1. 

Note that for x + and x_ > 1 one cannot clearly distinguish between 
first and second order terms. All of these limiting slopes agree with those 



given in || -||10||. However the F± in this paper are analytic functions with 
no unspecified constants. The functions 31F + and QF_ have been shown to 
continuously fit the first order dispersion data for the G-BN systems over 
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the whole range of values for x± J7|, and F± the dc conductivity results for 
Al x Gei- x 0. 

The slopes for QF + and for x± < 1 (second order terms) differ from 
those given in §| - which are based on the plausible assumptions made 
in |]J -||. However, no definitive experimental verification of these two terms 
seems to exist and the new experimental work given here strongly favours 
the second order term for 3?F_ given by Eqs.(|l]) and fl5|). The experimentally 
observed increases in 3?e with above C , some of which are given in Figs. 4 
and 5, is completely incompatible with the percolation equations, as given 
in - fli~0|l , but is qualitatively in agreement with Eq.(P if the separately 
measured dc parameters are used. Better agreement between Eq.(|I]) and 
experiment could be obtained if best fit parameters were used. 

From the above evidence Eq.([l|) may well be a better description of the 
ac and dc conductivity (dielectric constant) of percolative type systems than 
the standard percolation equations given in M - [IIJ. However, as Eq. (|T|) is 
a phenomenological equation its validity may need to be further tested by 
experiment as must the standard percolation equations. 
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Figure Captions. 
Fig. 1 

Plots of F + and F_ against x + {l)/uj+) and x-{oj/oj~), respectively. The 
parameters used are ~ C = 0.16, s = 1 and t = 2. As <fi ~ C the values 
used for a c must be accordingly large to yield a finite uf. The upper solid 
curve is $IF + and the lower one $IF_ (the second order dielectric loss term). 
The dashed line is ^sF_ (the first order term below C ) , and the dotted line is 
QF + . Note how this term rises above QF- in the region where x± is between 
1 and 100. 

Fig. 2 

A plot of the real part of the dielectric constant against frequency for a 
55%-45% G-BN powder on a log log scale for various values of (0 = 0.1309 
(open circles), 0.1290 (triangles), 0.1272 (open squares), 0.1236 (plus), 0.1219 
(crosses), 0.1203 (dots), 0.1187 (asteriks), 0.1171 (solid squares)). These are 
relative volume fractions, as the absolute error is about ±0.001. 

Fig. 3 
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A plot of the real part of the conductivity against frequency for a 55%- 
45% G-BN powder on a log log scale for various values of (0 = 0.1309 
(open circles), 0.1290 (triangles), 0.1272 (open squares), 0.1236 (plus),0.1219 
(crosses), 0.1203 (dots), 0.1187 (asteriks), 0.1171 (solid squares)). These are 
relative volume fractions as the absolute error is about ±0.001. 

Fig. 4 

Experimental values of 3?e for G-BN, plotted against 0, at lKhz (dots) 
and lMhz (circles). The nature and parameters for the theoretical curves are 
given in the text. 

Fig. 5 

Experimental values of 9?e for NbC coated grains, plotted against 0, at 
lKhz (dots) and lMhz (circles). The nature and parameters for the theoret- 
ical curves are given in the text. 
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